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Abstract 

Taking into account the additional operators coming up at finite temperature, we 
investigate the mass and residue of the nucleon in the framework of thermal QCD sum 
rules. We observe that the mass and residue of the nucleon are initially insusceptible 
to increase of temperature, however after a certain temperature, they start to fall 
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1 Introduction 


The investigation of the modifications on the hadronic parameters such as masses, decay 
constants, widths and coupling constants, etc. at hnite temperature is one of the impor¬ 
tant research areas of thermal QCD. The theoretical determinations of these modifications 
provide us with the better understanding of the properties of hot and dense QCD mat¬ 
ter produced by the heavy ion collision experiments. Such studies can also help us in 
understanding the internal structures of the dense astrophysical objects like neutron stars. 

In order to perform the theoretical studies on the hadronic parameters at hnite tem¬ 
perature, some non-perturbative approaches are required. One of the most applicable and 
powerful methods in this context is QCD sum rule approach. This method at zero temper¬ 
ature has been hrst proposed in [1, 2], which was based on the operator product expansion 
(OPE) and quark-hadron duality assumption. This approach has been extended to hnite 
temperatures in [3]. At hnite temperature, the Lorentz invariance is broken due to the resid¬ 
ual 0(3) symmetry. In order to restore the Lorentz invariance, the four-vector velocity of 
the medium is introduced. As a result, some new operators appear in the Wilson expansion 
[4-6] and the vacuum condensates are replaced by their thermal expectation values. 

In the literature, there are a lot of studies on the investigation of the hadronic pa¬ 
rameters at hnite temperature, but only a few of them have been devoted to the analysis 
of the modihcations of parameters such as the mass and residue of the nucleon at hnite 
temperature [7-14]. In [7], the nucleon mass and its residue as a function of the temper¬ 
ature have been investigated using the correlation function of a three-quark current with 
the help of the virial expansion. It has been obtained that the nucleon mass rapidly rises 
with temperature, while its residue falls with increasing temperature. In [8], the nucleon 
mass at hnite temperature has been analyzed via the QCD sum rules method. The authors 
have used lohe current and demonstrated that the nucleon mass is decrease with temper¬ 
ature. In [9], the nucleon mass shift at hnite temperature has been calculated using chiral 
counting arguments and a virial expansion. It has been shown that the nucleon’s mass 
decreases with increasing the temperature up to roughly T = 90 MeV, but after this point 
the mass starts to rapidly increase with increasing the temperature. In [10], the authors 
have obtained the spectral representation for the nucleon and given the expressions of the 
nucleon mass and its residue with respect to the temperature using lohe current by QCD 
sum rules approach. In this paper, S. Mallik and S. Sarkar have realized that the obtained 
results support the earlier results in [7]. In [11], the nucleon mass at hnite temperature has 
been calculated using Bethe-Salpeter equation (BSE) by observing that the nucleon mass 
decreases smoothly with temperature. In [12], the nucleon properties below the critical 
point temperature have been investigated. The author has solved the held equations in the 
mean-held approximation using the ehective mesonic potential at hnite temperature and 
found that the nucleon mass increases up to |Tc MeV then decreases near to the critical 
temperature T^. In [13], the nucleon mass at nonzero temperature is investigated within 
the framework of the linear a model. It has been shown that the nucleon mass increases 
monotonically with the temperature T. In [14], the nucleon properties at hnite temperature 
have been analyzed in the linear sigma model, in which the higher-order mesonic interac¬ 
tions up to the eighth-order interactions are included. It has been found that the nucleon 
mass increases with increasing of temperature. 
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In the present study, we calculate the mass and residue of the nucleon at hnite temper¬ 
ature using the most general form of the interpolating current in the framework of thermal 
QCD sum rules. In our calculations, we use the thermal light quark propagator in co¬ 
ordinate space containing the perturbative and non-perturbative contributions as well as 
additional operators rising at hnite temperature. We evaluate the two-point thermal cor¬ 
relation function hrst in coordinate space then transform the calculations to momentum 
space. In order to perform the numerical analysis, we hrstly determine the working regions 
of the auxiliary parameters entering the calculations and use the fermionic and gluonic part 
of the energy density as well as the temperature-dependent light quark and gluon conden¬ 
sates together with the temperature-dependent continuum threshold. We see that the Ioffe 
current which is a special case of the general current (t = — 1 ) remains roughly out of the 
reliable region. 

The outline of the paper is as follows: In the next section, taking into account the 
additional operators arising at hnite temperature, we derive thermal QCD sum rules for 
the mass and residue of the nucleon. The last section is devoted to the numerical analyses 
of the sum rules and comparison of the results with those existing in the literature. 

2 Thermal QCD sum rules for the mass and residue 
of the nucleon 

This section is dedicated to the details of the calculations of the mass and residue of the 
nucleon at hnite temperature using thermal QCD sum rules. The starting point is to 
consider the following two-point thermal correlation function: 


n(p,T) = i 


d*x e^P-^Tr(^p r(^J(x)J(0)^y 


( 1 ) 


where p is the four momentum of the nucleon, p = /Tre~^^ is the thermal density 

matrix of QCD at temperature T = l/P, H is the QCD Hamiltonian and T indicates the 
time ordering product. The most general form of the interpolating current for the nucleon 
can be written as 


J{x) 


2 r 


2 e, 


abc 


E 

2 = 1 


q'^’'"{x)CA\q^{x) 




( 2 ) 


where a, b, c are color indices, C is the charge conjugation operator. A] = I, Aj = A^ = 75 
and A 2 = t, with t is being an arbitrary auxiliary parameter and t = —1 corresponds to 
the lohe current (see [15-18] for detailed discussions on the nucleon current). Considering 
the isospin symmetry, we will deal with the proton (gi = u and q 2 = d) in the following. 

According to the general philosophy of the QCD sum rule approach, the above mentioned 
thermal correlation function can be calculated in two different ways. First, in terms of the 
hadronic parameters called the hadronic side and the second, in terms of QCD degrees 
of freedom called the OPE side. Equating the hadronic and OPE sides of the correlation 
function, the sum rules for the mass and residue of the nucleon at hnite temperature are 
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obtained. Applying a Borel transformation and continuum subtraction to both sides of the 
obtained sum rules, the contributions of the higher states and continuum are suppressed. 

In order to obtain the hadronic side of the correlation function, we insert a complete set 
of the nucleon state with the same quantum numbers as the interpolating current into the 
thermal correlation function. This side can be written in terms of different structures as 

n"“‘'(p,r) = nf“''(p",p„,r)^ + n"“‘(p",p„.r)^ 

+ n";‘‘(p2,p„,T)j(^ + n"“''(p2,p„,r)/, (3) 

where is the four-velocity of the heat bath, pq is the energy of the quasi-particle and I is 
the unit matrix. In the rest frame of the heat bath, = (1, 0, 0, 0) and = 1. To obtain 
the mass and residue of the nucleon, we use the structure whose coefficient is 

^"“*(p^p„,^) = -\Ut) ^ h,„, , (4) 

^ - m%{T) 

where niNiT) is the nucleon mass and \n(T) is its residue at hnite temperature. After 
performing the Borel transformation with respect to p^ in Eq. (4), the hnal form of the 
coefficient of on hadronic side is obtained as 

where is the Borel mass parameter. 

On the other hand, the OPE side of the thermal correlation function is calculated via 
OPE in deep Euclidean region. The correlation function on this side can also be written in 
terms of different structures as 

-h no^^(p2,po,T) + nf^®(p2,po,T)/. (6) 

To calculate the OPE side, we use the explicit form of the interpolating current in the 
thermal correlation function and contract out all quark pairs using the Wick’s theorem. As 
a result, we get 

n°'’®(p,r) = 

- (xhr.Tr Sf {x)ST'(x) '^+t{~„.Sf{x)^,ST'(x)S-f{x) 

+ Sf(p)Sr'(i)7,,Srh)75-75Sr'(i)rr Sf (x) 

- Sf(x)'„Tr Sf (p)Sf'(x)75 ) + (^Sf (xHS'p'(xh,Sf (x) 

- Sf(p)rr Sy(i)7,Sr'(i)7, )}) . (7) 
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where S' = CS'^C, Sq{x) is the thermal light quark propagator and Tr[...] stands for the 
trace of gamma matrices. The thermal light quark propagator is given as 
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1 — i—- 
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^ij 
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( 8 ) 


where rUq is the light quark mass, (gg) is the temperature-dependent quark condensate, 
is the external gluon held at hnite temperature, 0^^ is the fermionic part of the energy 
momentum tensor and are Gell-Mann matrices with A^A^^ = \{5ii5jk — -^SijSki). We 
insert the above quark propagator into Eq. (7) and perform the four-integral over x to 
go to the momentum space using Fourier transformations. To snppress the contribntion of 
the higher states and continnnm, a Borel transformation with respect to the as well as 
continuum subtraction are applied. We also use 


{Tr^Ga/sGfj,^) — —{ga^J.g|3u — gaug/3iJ.){Gx^G''^'^) 


1 r 


ga^gfiv gaugPfj, j^gPv ^a^ugUpU^gau T ^g^vgajjL) 


6 
X 


(9) 


to relate the two-gluon condensate to the gluonic part of the energy momentum tensor. Here 
is the traceless glnonic part of the stress-tensor of the QCD. After lengthy calculations, 
the Hp fnnction in Borel scheme np to linear terms in light qnark mass is obtained as 
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Finally, matching the hadronic and OPE sides of the correlation function, the sum rule for 
the residue of the nucleon at hnite temperature is obtained as: 

= 5nJ^^^(po,T). (11) 

To hnd the mass of nucleon from above sum rule, we apply derivative with respect to l/M^ 
to both sides of Eq. (11) and divide by itself. Eventually, the mass of nucleon at hnite 
temperature is obtained as: 


ml(T) 


-V) 


BUOP^{po,T) 


BUOP^{po,T) 


( 12 ) 


3 Numerical results 


In this section, we present numerical analysis on the mass and residue of the nucleon at 
hnite temperature obtained via thermal QCD sum rules. We determine how the results at 
hnite temperature shift from those obtained at zero temperature. For this aim, we use the 
input parameters given in Table 1. 

Besides these input parameters, we need to know the temperature dependence of the 
quark and gluon condensates, the thermal average of the energy density and the temperature- 
dependent expression of continuum threshold. 

In this study, for the temperature-dependent quark condensate we use the parametriza- 
tion given in [23], which can be htted to the following function up to the quark-gluon 
deconhnement at a critical temperature Tc ~ 197 MeV: 


{qq) 


_ {o\qq\o) _ 

X -)- gl8.10042(1.84692r2+4.99216r-l) ’ 


(13) 
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Table 1: The values of input parameters used in numerical calculations [19-22], 


where (0|gg|0) denotes the light-quark condensate in vacuum. The above parametrization 
is consistent with the lattice QCD results [24,25], 

For the gluonic and fermionic parts of the energy density we use the results obtained 
from lattice QCD at higher temperatures. In the rest frame of the heat bath, the results 
of some quantities obtained using lattice QCD in [26] are well htted by the help of the 
following parametrization: 

(05p) = (0^p) = 2-4g(ii3.867r2-i2.i90r) _ 10.141T^ (14) 


where this is valid for T > 130 MeV. 

We take the temperature-dependent continuum threshold as 


so(T) ~ So 


{m) 

(0|gg|0)’ 


(15) 


where Sq is the continuum threshold at zero temperature. The temperature-dependent 
gluon condensate corresponding to Tc ~ 197 MeV is also given as [27] 


(G^) = (OIG^IO) 


1.65(^)®-^^® + 0.04967(^)°-^2^^ 

-L r J- r ■ 


( 16 ) 


where (0|G^|0) stands for the gluon condensate in vacuum. 

The next step is to determine the working regions for three auxiliary parameters, namely 
Borel mass parameter M^, the continuum threshold sq and mixing parameter t. It is 
expected that the physical quantities are roughly independent of these auxiliary parameters 
in their working regions. 

The working region for the Borel mass parameter is found demanding that the con¬ 
tributions of the higher states and continuum are adequately suppressed as well as the 
contributions of the higher dimensional operators are small. These conditions lead to the 
interval 0.8 GeV^ ^ ^ 1.2 GeV^ for the Borel mass parameter. The continuum thresh¬ 

old is not utterly arbitrary but it is associated with the energy of the hrst excited state 
with the same quantum numbers as the chosen interpolating current. From our numerical 
analysis, we obtain the interval sq = (1.5 —2.0) GeV^ for the continuum threshold. In order 
to determine the working region of the mixing parameter t, we need to investigate the de¬ 
pendence of the physical quantities on t in the whole range —cxd < t < oo. For this purpose, 
we plot the mass and residue of the nucleon with respect to a; = cos6, where t = tanO, 
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at fixed values of the continuum threshold sq for Borel mass parameter = 1 GeV^ at 
T = 0 in hgure 1. From this hgure, we see that the dependences of the mass and residue 
on this parameter are relatively weak in the intervals —0.6 < x < —0.2 and 0.2 < x < 0.6. 
Also, we realize that the Ioffe current which corresponds to a: = —0.71 remains out of the 
reliable region. 




Figure 1: The mass and residue of the nucleon as a function of x at = 1 GeV^, T = 0 
and at different hxed values of sq. 

To see how our results depend on the Borel mass parameter in the above determined 
working regions, we plot the mass and residue of the nucleon at zero temperature with 
respect to Borel mass parameter at different hxed values of x and sq in figure 2. This 
figure depicts that the values of the physical quantities under consideration depend weakly 
on the Borel parameters in its working region, 0.8 GeV^ ^ ^ 1.2 GeV^. 



MfiGeVfi M^[GeVh 


Figure 2: The mass and residue of the nucleon as a function of at T = 0 and at different 
fixed values of Sq and x. 

In further analysis, the dependence of the mass and residue of the nucleon on the 
temperature are shown in figure 3 at different values of the auxiliary parameters in their 
working regions. With a quick glance at this figure, we observe that the mass and residue 
of the nucleon remain approximately unchanged up to T = 0.15 GeV, after which they 
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Figure 3: The mass and residue of the nucleon as a function of temperature at = 1 GeV^ 
and at different hxed values of Sq and x. 
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Figure 4: Mn/Mn{0) and Xn/^n{0) as a function of T/Tc at = 1 GeV and at different 
hxed values of sq and x. 

start to diminish increasing the temperature. Near to the critical temperature, the mass of 
the nucleon falls to roughly 42% of its values at zero temperature while the residue reaches 
roughly to 60% of its value at zero temperature. 

To better see the deviations of the mass and residue of the nucleon from their vacuum 
values, we show the dependence of the ratios M]\f/M]y{0) and AAr/A 7 v( 0 ) to T/Tc in hgure 
4. When we compare these behaviors in terms of the temperature with those of existing 
results in the literature, we see that as far as the residue is concerned our results support 
the behaviors obtained in [7,10] for the residue with respect to the temperature at whole 
region. In the case of mass, our results on the variation of the mass of the nucleon are in 
agreement with those of [8,11]. Although are also in good agreements at low temperatures, 
our results on the behavior of the mass with respect to temperature are in contradiction 
with those of [7,9,10,12-14] at high temperatures. 

To summarize, we calculated the mass and residue of the nucleon at finite temperature 
in the framework of the two-point thermal QCD sum rules. We used the general form of 
the nucleon’s interpolating current and found the working regions of auxiliary parameters 
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entering the calculations. It has been obtained that the Ioffe current remains out of the re¬ 
liable region for the general mixing parameter entering the general form of the interpolating 
current. We have also used the temperature dependence of the quark and gluon conden¬ 
sates, the thermal average of the energy density and the temperature-dependent expression 
of continuum threshold to numerically analyze the sum rules for the mass and residue of 
the nucleon. We found that the mass and residue of the nucleon remain unchanged up to 
roughly 0.15 GeV. After this point, they start to decrease with increasing the tempera¬ 
ture such that near to the critical temperature they fall down, considerably. This can be 
consider as a signal for transition to the quark-gluon plasma phase. We also compared the 
behavior of the mass and residue of the nucleon in terms of temperature obtained in the 
present work with those existing in the literature. 
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